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The possible quantum Hall ferromagnet at a filling factor v = is investigated for 
the zero-energy (N = 0) Landau level of the two dimensional massless Dirac fermions 
in a-(BEDT-TTF)2i3 under pressure with tilted cones and a twofold valley degener- 
acy resulting from time-reversal symmetry. In the case of the Dirac cones without 
tilting, the long-range Coulomb interaction in the N = Landau level exhibits the 
SU(2) valley-pseudo-spin symmetry even to the order 0(a/lu), in contrast to N ^ 
Landau levels, where a and /h represent the lattice constant and the magnetic length, 
respectively. Such a characteristic comes from a fact that zero-energy states in a par- 
ticular valley are restricted to only one of the spinor components, whereas the other 
spinor component is necessarily zero. In the case of the tilted Dirac cones as found in 
a-(BEDT-TTF)2l3, one obtains a non-zero value of the second component and then 
the backscattering processes between valleys becomes non-zero. It is shown that this 
fact can lead to easy-plane pseudospin ferromagnetism (XY-type). In this case, the 
phase fluctuations of the order parameters can be described by the XY model leading 
to Kosterlitz-Thouless transition at lower temperature. In view of these theoretical 
results, experimental findings in resistivity of a-(BEDT-TTF)2l3 are discussed. 

KEYWORDS: massless Dirac fermion, Weyl equation, zero-gap state, a-(BEDT-TTF) 2 l3, graphene, 
magnetic field, Landau level, electron correlation, quantum Hall ferromagnet, Kosterlitz- 
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§1. Introduction 

The massless Dirac fermions in the quasi- two-dimensional organic conductor a-(BEDT- 
TTF) 2 l3, 1 ^ which obey the tilted Weyl equation 2 ' 3 \ have attracted much interest because 
of the mysteries of the experimental findings such as the weak temperature (T) dependence 
of resistivity (close to h/e 2 ), the strong T-dependence of the Hall coefficient, 4 " 7 '' and the 
two-step increase of resistivity with decreasing T in the presence of the magnetic field. 5 -* 

In the absence of magnetic field and under pressure P = 18kbar, the in-plane resistivity 
decreases weakly with decreasing T from the room temperature and turns to increase below 
the onset temperature, 10K, and then it saturates in the limit of T — > O. 5 -* The origin 
of such weak but obvious T-dependence has not been elucidated yet. In the presence of 
magnetic field, H, perpendicular to the conducting plane, the two-step increase of resistivity 
is observed. 5 '' For example at H = 10T, with decreasing T, the resistivity decreases weakly 
from room temperature and turns to increase below T = 20K with the plateau in the lower 
temperature region. Another relative sharp increase is observed around T} = 5K leading to 
apparent saturation as T — ► 0. Both T and T} increase with increasing magnetic field. 

In the absence of tiling it has been established that the energy spectrum of the massless 
Dirac fermions becomes discrete by the Landau quantization owing to the orbital motion in 
the magnetic field, = sgn(N) J 2%v 2 eH\N\/c, where v is the velocity, and N is an integer. 
Each Landau level has large degeneracy proportional to H, and is split into two states with 
up and down spins by the Zeeman energy, as shown in Fig. 1. The effects of tilting have 
been investigated recently and it is found that the Landau-level structure is qualitatively the 
same and that the effects of tilting are incorporated as modifications of effective velocity. 
By use of values of relevant parameters appropriate for a-(BEDT-TTF)2i3, 3 ' 8 ' ) E\ = 5meV 
and 2Ez = g^H = lmeV with the g-factor g = 2 for H=10T. Hence we see that the energy 
scale seen in resistivity measurement, To (= 2meV) and T, are smaller than E\. Thus, 
the observed two-step increase of resistivity in a-(BEDT-TTF)2i3 may be attributed to the 
N = Landau levels, whose causes will be studied theoretically in this paper. 

The long range Coulomb interaction plays an important role for massless Dirac fermions. 
The effective Coulomb interaction under magnetic field, /, is estimated as / = e 2 /e/n — 
50yjH [T]/e meV, where In is the magnetic length, l# = ^Jfic/eH. Although the polarizability 
e of a-(BEDT-TTF)2l3 under high pressure has not been identified so far and then there is 
some ambiguity, it is demonstrated later that = Landau states have instability toward 
the pseudo-spin ferromagnetism since the effective Coulomb interaction I can exceed 2E Z . In 
the presence of tilting, it is shown that the electron correlation can give rise to the quantum 
Hall ferromagnet of the pseudo-spin (the degree of freedom on the valleys) with the help 
of the large degeneracy of the Landau levels. The easy plane anisotropy of the pseudo-spin 
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ferromagnet results from the back scattering processes which is the inter-valley scattering 
terms exchanging large momentum. Moreover, it is shown that the effects of fluctuations of 
phase variables of the order parameters can be described by XY Heisenberg model leading 
to Kosterlitz-Thouless (KT) transition 9 - 1 at lower temperature. 




Fig. 1. Schematic figure of the Landau levels as a function of magnetic field. The energy scale of the onset 
temperature of the first increase of resistivity, To, is located between the Zeeman gap 2Ez and the E\ — Eq 
gap, where the Boltzmann factor fcs is taken as unit. The energy scale of the second increase of resistivity 
at low temperatures, T, is smaller than the Zeeman gap. 



§2. Formulation 

2. 1 Hamiltonian for massless Dirac fermions with tilting 

In the absence of magnetic field, the Hamiltonian of the massless Dirac fermions is given 
by, 

Ti = + 

k77'oT 

H' = -J J drdrV (r - r>(r)n(r') (1) 

with the long-range Coulomb interaction Vo(r) = e 2 /er and the density operator n(r). 
The degree of freedom on the spins are represented as a = ± corresponding to f , I- The 
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degree of freedom of the pseudo-spins, r = ±, corresponds to the valleys R, L. The valleys 
are located at the crossing points of the conduction and valence bands, ±k , where k is 
an incommensurate momentum in the first Brillouin zone. 1 ) The creation and annihilation 
operators, c\. and Ck 7crT , respectively, are based on the Luttinger-Kohn representation 10 ^ 
using the Bloch's functions at the crossing points as the basis of wave functions, 2 ) and then 
7 = 1,2 denotes the basis of the Luttinger-Kohn representation. The relation between the 
Luttinger-Kohn representation and the site representation based on the molecular orbitals 
are described in Appendix. The low-energy properties around the two crossing points, labeled 
by r = ±, are described in terms of the two tilted Weyl Hamiltonians 2 ' 3 -* 

Hq' t ~ + = h(v~k ■ <t + w ■ ko"o) 

H°' T= ~ = -H(vk ■ (T* + w ■ ka ) (2) 

with respect to the time-reversal symmetry, where v represents the velocity of the cone and 
wo represents the tilting velocity. Here, we have neglected the anisotropy of the velocity 
of the cone. In the case of the cone with anisotropic the velocity, one needs to replace 
uk > (v x k x ,Vyky). 

Once a magnetic field perpendicular to the conducting plane is taken into account, the 
momentum k is replaced by k + (e/c)A with the vector potential in the Landau gauge 
A = (0, Hx, 0), and then we obtain 



^0 



TT<T,T = 

n 



Hw 0x \£ + w 0y {\l + oE z hv{\-§- x - itf « + £JE)} 

19 i „„ (Id | eHx\\ „rp *„,r 19 \(1 d , eHx 



Hw 0x \f x + w 0y (\% + - aE z hv{-\& - + 

in terms of the Zeeman energy E z . 
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2.2 Zero-energy Landau level for the case with tilting 
The eigen equations for the zero-energy Landau level are given by 

HS^ t x (t)=0 (4) 

which gives the eigen functions of the tilted Weyl Hamiltonians in the presence of magnetic 
field. The wave functions are given by 3 ) 

#e(r) = -^=e- iXy/l H V T {x - X)e- iTk(rr (5) 

with 
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fix) = . X Xre^ x2/21 -, (6) 
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where X is the guiding-center coordinate and L is the length of the system. The spinor parts 
Xt ar e given by 

XT = + /*o 2 + (l + 7) 2 V 

Xr--= , 1 ( | • (8) 

Iwl + il + ^f \ -woe** 



Here, we have defined 



= ^±^y, ( 9 ) 



in terms of the effective tilting parameter 



w 0x \ 2 . fw 0y x 



™o = \l(^)+{^) (10) 



with w = (w 0x , w 0y ), and 7 = y 1 — iDq. We note that one recovers the usual result for the 
n — wave function in graphene when we use in the limit wq — > (7 — > 1), i.e. in the case 
without tilting. 

2.3 Effective Hamiltonian on fictitious magnetic lattice 

We consider the present model by using the bases of the Wannier functions for the magnetic 
rectangular lattice which is introduced fictitiously. The wave functions in the Landau gauge 
are not localized in the y-direction, but in the x-direction around the position X. By 
applying the periodic boundary condition in the ^/-direction with the system length L, the 
X is discretized as X = —2Txl\jjL with an integer j. The Wannier functions, which satisfy 
orthonormality and are localized around R = (ma, nb) with integers m and n as shown in 
Fig. 2, are constructed by the linear combination of <^x( r )> 

*r(*) = ^1 (*'] dX e iXnb / l «P ( r) j (ii) 
a\/b J-a/2 

where a is arbitrary but a ^> 2i\l\jL and b = 2i\l\la. We note that |$^(r)| exhibits 
an exponential-like decrease in the x-direction, but decreases algebraically as \y\~ l in the 
y-direction. 

In the basis of these Wannier functions, the effective interaction Hamiltonian is given by 

n' N=0 = l EE/ /drdrVo(r-r') 

R-1,2,3,4 CO 1 ' 71,2,3,4 

x ^(^-^(r)] [^3(0^^(0; 

XC R.lcrri C R20-T2 C R3 -' r3 C R40-'r4 (12) 
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Fig. 2. The r-dependence of |$^(r)| 2 with R = 0, where $R_(r) is the Wannier functions for the magnetic 
rectangular lattice. 



under the assumption that the density operator is effectively determined by the field operator 
for N = Landau states and the contributions to the Hamiltonian from N ^ states are 
negligible, i. e. we use the density operator 

n(r)= £ <^)]cl aT c WaT , (13) 

RR'crrr' 

Thus the effective Hamiltonian for the N = Landau states in the magnetic rectangular 
lattice is given by 

ftcff = J2(-crEz)cl aT c iaT 

1.(7 T 

+ ^2 ^ijkl C laT C jo-rc\. a / T fCl a ' T ' 
ijklacr'TT 1 

+ Yl W ijklcLfCjaTcl a , T C la > f , (14) 
ijklaa'r 

where i, j, k, and I denote the unit cells of the magnetic rectangular lattice at R i; R J; R&, 
and R/, respectively, and f = — r. 
The forward-scattering term, Vijki, is given by 
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V m =t.l /drdr'Vo(r-r') $T( r )t . $J( r ) <( r ') f • $['(r') (15) 



2 

from the long wave length part o{7i' N=0 . This term does not depend on the spin and pseudo- 
spin, and then it does not break the SU(4) symmetry, neither in the spin subspace nor in 
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that of the pseudo-spin. We find that the forward-scattering term is not affected by the 
tilting, because {x\ ■ Xt){x\> • Xr') = 1- 

On the other hand, the backscattering term, Wijki, which is the inter- valley scattering 
term exchanging large momentum 2ko and breaks the SU(2) symmetry in the subspace of 
the pseudo-spin, is given by 

W ijkl = \j j drdr'V (r - r') [$[(r)t • $J(r)] [^(r') 1 • $[(r')] (16) 

from the short wave length part of 7~L' N=0 . In the absence of tilting, as e.g. in graphene, the 
backscattering term vanishes because (x\ • Xr) = 0. 12 ** We find that the tilting is essential to 
have a non-zero backscattering term. The tilting dependence of the backscattering term is 
given by the spinor part, 

(xt ■ Xr)(xt ■ X.) = { ^J\ 1+ ^ )2]2 - + 0(w*), (17) 

where the last step has been obtained from the limit Wo <C 1. The ratio between the forward 
and the backscattering terms, Wijki/Vyki, is given by 



Wijki/Vijki 



In 7 

where ol is the lattice constant in the conducting plane. The backscattering term is pro- 
portional to ol, since the large momentum |2/c | — tt/ol is exchanged. 12 -* We note that the 
lattice constant of a-(BEDT-TTF) 2 l3, ol — 10 A, is much larger than that of graphene. Thus 
it is expected that the backscattering term plays an important role for electron-correlation 
effects in a-(BEDT-TTF)2i3. The typical value of the ratio Wijki/Vijki is approximately 0.07 
for a-(BEDT-TTF)2l3 at H = 10T using the tilting parameter t^o — 0.8. The Umklapp scat- 
tering term (rfrf -term) can be neglected, because it is exponentially smaller than the other 
terms as a function of ol/^Hj which is estimated as 0.1 at H = 10T in a-(BEDT-TTF) 2 l3. 

§3. Pseudo-spin ferromagnet and Kosterlitz-Thouless transition 



Possible spin and pseudo-spin ferromagnetic states in the zero-energy Landau level in 
graphene have been extensively studied in recent years. 12_19 ) Generically, the ferromagnetic 
ordering may be understood within an interaction model with no explicit spin or pseudo-spin 
symmetry breaking; in order to minimize their exchange energy, the global iV-particle wave 
function should be fully antisymmetric in its orbital part, the (pseudo-)spin part needs to 
be fully symmetric in order to fulfil fermionic statistics. Whereas in a normal metal this 
ordering is only partial, due to the increase in the kinetic energy, a single Landau level may 
be viewed as an infinitely flat energy band, and the ferromagnetic ordering may therefore be 
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complete. In the absence of an explicit symmetry breaking, such as the Zeeman effect that 
naturally tends to polarize the physical spin or the above-mentioned backscattering term 
that affects the pseudo-spin, no particular spin or pseudo-spin channel is selected, and one 
may even find an entangled spin-pseudo-spin ferromagnetic state. 20 ) The symmetry-breaking 
terms may, thus, be viewed as ones that choose a particular channel (spin or pseudo-spin) 
and direction of a pre-existing ferromagnetic state by explicitly breaking the original SU(4) 
symmetry. 



3.1 Mean-field solution 
The mean-field Hamiltonian for the pseudo-spin ferromagnetic state is given by 



n 



MF 



E 



j err 



-ge z -X!* r 'V./'""'. + Wt 



C jar C jo-r 



"2 E ( V W + W Hjj)(( C LTCiaf)c] af C jr7T + h.c) 



(19) 



l]OT 



and the order parameter of the pseudo-spin ferromagnetic state, A, which is independent of 
% and a, is defined by 

A = 2/<4_< W ) (20) 

with the effective interaction / = J2i(Vmi + Wuoo). Using the spin polarization, m = 
J2ar °~ n ioTi which is also independent of i, and the renormalized Zeeman energy, Ez = 
Ez + ml, the mean-field Hamiltonian is given by 



H 



-E z 


-A* 








-A 


—E z 














E z 


-A* 








-A 


Ez 



with c~ 



[CjT+> Cjf_, Cji + , Cjl- 



). The mean-field solution is calculated from 



21 r 

A| = — / dx/(x)Im- 

7T J 



IAI 



x + i5 - aE z ) 2 



IAI 



and 



m 



J2af(-aE z + r\A\) 



(21) 



(22) 



with the Fermi distribution function, f(x). 

The ground state in the case with / < Ez is a spin polarized state without pseudo-spin 
polarization (m = 1 and A = at T = 0), where electrons reside in the spin-up branches of 
the N = Landau levels (see the left hand side of fig. 3(a)). The ground state in the case 
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with I > Ez is, on the other hand, a pseudo-spin ferromagnetic state (m = and |A| = I 
at T = 0), where the easy-plane pseudo-spin polarization lifts the pseudo-spin degeneracy 
and then the spin polarization is suppressed (see the right hand side of fig. 3(a)). 

The mean-field phase diagram in the I-T plane scaled by Ez is shown in fig. 3(b). The 
transition temperature for the easy-plane pseudo-spin ferromagnetic state, T c , is finite in the 
case with I > Ez, and increases with increasing /. Below T c , the spin polarization vanishes 
at T — > 0, although it is still finite at finite temperatures. 

One notices that this competition between a spin polarized state and an easy-plane pseudo- 
spin ferromagnetism is original to the filling factor v = 0, where necessarily two (of four) 
subbranches of the zero-energy Landau levels are occupied. In contrast to this particular 
filling factor, this competition is absent at v = ±1, where only one subbranch is occupied and 
where, therefore, a spin polarization does not exclude a simultaneous pseudo-spin ordering 
in a coherent superposition of both pseudo-spin states. 



E z +ml+|A| 




Fig. 3. (a) Schematic figure of the energy levels in the spin-polarized state (left hand side) and the pseudo- 
spin ferromagnetic state (right hand side), (b) Phase diagram in the plane of the interaction / and 
temperature T scaled by the Zeeman energy Ez- 



3.2 Phase fluctuations and Kosterlitz-Thouless transition 

In the presence of an order parameter with a finite amplitude below T c , phase fluctuation 
exists with the characteristic length of spatial variation much longer than the fictitious lattice 
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spacing. The effect of these phase fluctuation, which has so far been ignored in the mean- 
field approximation, is treated on the basis of the Wannier functions and the resulting model 
is similar to the XY model leading to the KT transition. Using the pseudo-spin operator, 
§£_ = c\ aL Ci a R, and the real spin operator, Sj Z = ~ yj CTT aCj aT Cj aT , the mean-field Hamiltonian 
is given by 

H M F = -2(E z + mI)J2S jz 

j 

- 2 J2 hi ((Sl)SJ + + (S° + )SJ_) . (23) 

ijcr 

The real spin polarization, m, remains finite at finite temperatures, although it vanishes 
at T — > in the pseudo-spin ferromagnetic state. However, it can be shown that the spin 
polarization is independent of i, since the pseudo-spin can fluctuate only in the easy plane 
and then the occupation numbers of electrons are independent of i. The interactions between 
the pseudo-spins on the magnetic rectangular lattice, 7y = + Wajj. The interaction 1^ 
rapidly decreases with increasing |Rj — Rj|. It is numerically found as seen in Fig. 4 that 
the nearest-neighbor 1^ is approximately isotropic and the ratio I^i+i/I^i = 0.10 with the 
arbitrary choice of b = y/2a (but leading to the almost isotropic localization of the Wannier 
function on the fictitious lattice), where 

'iiii ^^iiii 'iiii 

The phase of the order parameter corresponds to the angles of the pseudo-spin. The x- 
and ^-components of the pseudo-spins are given by (S^) = Re(Sf_) and (S 1 ^) = — lm(S?_), 
respectively, with (S^ + ) = (Sf_)*. Thus (S£_) can be represented using the angle of the 
pseudo-spin from the x-direction in the x-y plane, as (Sf_) = | (S"_) | exp(— i^f). When 
the characteristic length of spatial variation of the phases is much longer than the lattice 
spacing, we can expand the free energy by the fluctuations of the phases — 1 — cos(0J — 
4>1) <C 1 under the assumption that the amplitude, \(S^_)\, does not change within the 
characteristic length of the phase fluctuation, and then we obtain 

Ff -F n ^(F f -F n )\ f ^ 0+ £ d{Ff - Fn) f% 

= F - £ ^cos(0J-C) (25) 

with J[j = A\(SZ.)\ 2 Iij, where Ff and F n denote the free energies of the pseudo-spin ferro- 
magnetic and normal states, respectively, and Fq is independent of the phases. The effects of 
the phase fluctuations in Jf- are neglected since those effects are the higher-order terms, and 
J?j — > Iij at temperatures much lower than T c . Then the physics of the phase fluctuation is 



10 



equivalent to that of the two-dimensional XY Heisenberg model with the nearest-neighbor 
exchange interaction, J = i^+i = 0.087/. In the two-dimensional XY model, the KT tran- 
sition occurs due to the onset of bound pairs of the vortices. 9 ^ It is known that Tkt — 1-54 J 
by the renormalization group analysis. 21 ' Here the effects of spin polarization are negligible, 
because J is much larger than the Zeeman energy and the phase fluctuation on each electron 
with up or down spin is described by the same interaction, J. 

Lastly we discuss the role of the long-range part of 1^ farther than the nearest neighbor 
one. Figure 4 shows the distance dependences of I 0i j/I 0;0 with Rj = (na, 0) (the closed 
circles) and Rj = (0, nb) (the open circles) defined on an integer n, where we take b = \[2a. 
The interaction Iqj decays very rapidly along the x-axis but slowly along y-axis. Then 
the role of the long range part of Jjj along y-axis should be considered for improving the 
effective Hamiltonian. The KT transition temperature Tkt, however, is determined by the 
competition between the excitation energy of a vortex and the entropy effect coming from 
the degree of freedom for the position of the vortex core. Since the length scale of the vortex 
is much longer than a and the interaction is ferromagnetic, the long-range part of Jy does 
not disturb the KT transition essentially. 



10 : < c 
10 ? 

o 
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Fig. 4. The n-dependcnccs of Io,j/Io,o with Rj = (na, 0) (the closed circles) and Ri = (0,nb) (the open 
circles), where we take b = \/2a. The dashed and dotted lines are a guide to the eye. 
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§4. Relation between Experimental Findings and Theoretical Results 

In the presence of magnetic field perpendicular to the conducting plane, the two-step 
increase of resistivity, for example at T = 20K and T} = 5K at H = 10T, is observed in 
a-(BEDT-TTF)2l3, where both To and T\ increase with increasing magnetic field. 5 ) We may 
be able to associate two stepwise changes of resistivity as due to the easy-plane pseudo-spin 
ferromagnetic transition at T c and the KT transition at Tkt, since our estimate indicates 
that T c = 4Tkt, where T c = 0.5/ as seen in Fig. 3(b) in the region of I/Ez » 1 of interest 
and Tkt — 1-54J = 0.137. We emphasize that the tilting of the Dirac cone is essential to 
the appearance of the easy-plane pseudo-spin ferromagnet, and thus, to the appearance of 
the KT transition, due to the long range Coulomb interaction. 

§5. Conclusion and Discussion 

In the present paper, motivated by the experimental observation of the particular temper- 
ature dependences of resistivity in a-(BEDT-TTF) 2 i3 under magnetic field, the possibility 
of the pseudo-spin quantum Hall ferromagnet at v = has been investigated in the massless 
Dirac fermion system. The pseudo-spin ferromagnetic transition occurs when the electron 
correlation exceeds the Zeeman energy. The tilting of the Dirac cone induces the backscat- 
tering terms resulting in the easy-plane pseudo-spin ferromagnet. There will be intrinsic 
fluctuations and T c should be considered only as a crossover temperature for the growing am- 
plitude of order parameters with remaining large phase fluctuations in the two-dimensional 
system. To treat such phase fluctuations, a spatially localized basis set similar to 'Wannier 
function" are introduced, which indicates that the model is similar to the XY model which 
is known to lead to the KT transition at lower temperature, Tkt- I n comparison with ex- 
periments, the two-step increase of resistivity with decreasing temperature are observed at 
around 20K and 5K at H = 10T. Present theory has revealed Tkt — 5K on the choice of the 
parameters giving T c = 20K, and then there are reasonable correspondences to identify two 
stepwise changes of resistivity as due to the amplitude growing and the phase coherence of 
the order parameters. 

Obviously there are remaining problems to be clarified. The experimental data indicates 
the saturation of resistivity in the low temperature. The saturation indicates the existence 
of dilute carriers which may originate from weak three-dimensionality or disorder. 

In graphene, 22 ** the the quantum Hall ferromagnet at v — has been investigated, 
12,14,13,16,19) anc j ver y reC ently it is suggested that the electron-phonon interaction break- 
ing the pseudo-spin SU(2) symmetry, which may be characteristic of graphene, induces the 
easy-plane pseudo-spin ferromagnet resulting in the KT transition 23 ^ in order to explain the 
possible KT transition observed in graphene. 24 ** We emphasize that the backscattering term, 
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which is the key factor for the easy-plane pseudo-spin ferromagnet in our paper, is charac- 
teristic of a-(BEDT-TTF) 2 l3, but can be realized in graphene by distorting the honeycomb 
lattice. In addition, we note that a lattice model describing the fluctuations of the pseudo- 
spins should be based on Wannier functions which satisfy orthonormality under magnetic 
field, since the bases on the original crystal lattice are no longer the eigenstates of the N = 
Landau levels. 

The effects of the short range parts of the Coulomb interaction in graphene also have 
been investigated. If once the pseudo-spin ferromagnet ism occurs, the Hubbard-£/-type on- 
site interaction favors the easy-plane ferromagnetism (uniform charge density), while the 
nearest-neighbor interaction V favors the easy-axis ferromagnetism (the sublatice CDW), 
in the case of graphene. 14 -* In a-(BEDT-TTF)2i3, however, the easy-axis ferromagnetism 
does not correspond to CDW directly, because the bases of the Weyl Hamiltonian is not the 
sublattice but the Bloch states at k = ±k . The Bloch states at k = ±k consist of the 
linear combination of the contributions from four BEDT-TTF molecules. Thus, although the 
easy-axis ferromagnetism may modify the intrinsic charge disproportionation in a-(BEDT- 
TTF) 2 l3, the effect of V on such state may be weaker than than of graphene. It is an 
interesting difference between graphene and a- (BEDT-TTF) 2I3, and it will be investigated 
intensively in future. 

Lastly, we discuss the renormalization of fluctuation in the pseudo-spin ferromagnet, which 
are very complicated and not captured on the mean-field level. In the absence of the Zeeman 
effect, the ferromagnetic moment may fluctuate in the SU(4) space at temperatures between 
T c and a symmetry breaking temperature, T s b, which is essentially given by the symmetry- 
breaking interaction energy, W^ki- The ferromagnetic moment may be forced in the easy- 
plane below T s b- In the presence of the Zeeman effect, the situation is more complicated 
owing to the competition between Wijki and the Zeeman energy, Ez- The results in the 
present paper, thus, identify a new research target, i. e. a two-dimensional SU(4) model 
with the symmetry-breaking terms. 
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Appendix: Relation between Luttinger-Kohn Representation and Site Repre- 
sentation 

A basic model describing the two-dimensional electronic system in a-(BEDT-TTF)2i3 is 
given by 25 " 27 ) 

W s ite = (tiorjP a taa a jP<r + n - C -)) (A' 1 ) 

(ia:jf3),cr 



11 



where i and j denote the unit cells, a ,(3 (=A, A', B and C) denote the molecular orbital 
sites in a unit cell, a\ aa denote the creation operators on the site representation, and U a .jp 
is the transfer energy between (i,a) site and (j,(3) site. Using the Fourier transformation, 
we obtain 

H-site = X/ e a/3cr(kO a kaGr a k/3o- 
ka/3<r 

e af3 (k) = Y / t a( 3^ 5 , (A2) 

5 

where 5 denotes the vector representing the nearest neighbor of the unit cell. Such Hamil- 
tonian (here we call it the site Hamiltonian) is diagonalized by the eigenvalue equation 

4 

e a/ 3 CT (k) dp^ik) = £ 7CT (k) d a ^(k), (A-3) 

(3=1 

where ^ 7<T are the eigenvalue with the descending order, £io-(k) > £20- (k) > £30- (k) > ^ 4cr (k), 
and d ara (k) (7 = 1,2,3,4) are the corresponding eigenvectors. Here £io-(k) and ^2<r(k) are 
the conduction and valence bands, respectively, because there are six electrons in the four 
molecules, i.e., the 3/4- filled electronic system. 

Expanding the site Hamiltonian in the linear order of momenta from ±k[, and using the 
Luttinger-Kohn representation based on the Bloch's functions at ±k[, (corresponding to 
t = ±), we obtain the tilted Weyl Hamiltonians, where ±k[, are infinitesimally close to the 
crossing points ±ko, respectively. 2 ' 28 -* Then the creation operators on the Luttinger-Kohn 
representation cj C7crr are given by 

cUar = H < Cy<T( rk o) Lr- ( A ' 4 ) 
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